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ABSTRACT 

Flow structure of self-excited. laminar, axisymmetric, 
momentum-dominated helium jets discharged vertically 
into ambient air was investigated using high-speed 
rainbow sehlieren defleetometiy technique. 
Measurements were obtained at temporal resolution of 
I ms and spatial resolution of 0.19 mm for two test 
cases with Richardson number of 0.034 and 0.0 1 8. 
Power spectra revealed that the oscillation frequency 
was independent of spatial coordinates, suggesting 
global oscillations in the (low. Abel inversion algorithm 
was used to reconstruct the concentration field of 
helium. Instantaneous concentration contours revealed 
changes in the flow field and evolution of vortical 
structures during an oscillation cycle. Temporal 
evolution plots of Itclium concentration at different 
axial locations provided detailed information about the 
instability in the flow field. 

Keywords: Color Sehlieren. Momentum-dominated 
jets. Flow Instability. 

INTRODUCTION 

Flow structure and instability of gas jets has been the 
subject of several investigations due to their wide 
ranging applications. When a low density jet is injected 
into a high density environment, highly periodic 
oscillations occur in the flow field under some 
conditions. Tlie important operating parametets in the 
previous studies included the jet Reynolds number 
(Rc-d.UiV). the jet Richardson number iRi-g.d,(p.- 
pVlp.U')). the momentum thickness (0) and the density 
ratio (S=(Vp,|. The instability has been described using 
flow -oscillation frequency (0 nr the Strouhal number 
(St-f.dU). Low -density jets can be classified as cither 
buoyancy dominated (Ri-T> or momentum dominated 
<Ri«l). In the literature, we find several studies 
conducted in both jet flow regimes. 
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Subbarao and Cantwell' identified the oscillating and 
non-oscillating regimes for vertical liclium jets in the 
Reynolds number-Richardson number space using a 
stroboscopic sehlieren technique. Measurements 
revealed large centerline velocity fluctuations, and early 
and abrupt breakdown of the potential cote. They 
reported the oscillating behavior of the helium jet at 
moderate v alues of Richardson numbers (0.5 < Ri < 6) 
and stated that this type of flow was subjected to an 
unusual type of transition to turbulence consisting of a 
rapid but highly stmetured and repeatable breakdown 
and intermingling of the jet with live free stream fluid. 
The strong dependence of Strouhal number on 
Richardson number indicated the dominance of 
buoyancy efl'ect. 

Mamins et al. : studied oscillations in both buoyant 
helium jets and flames. They used shadowgraph 
technique to observe the oscillation frequencies of 
helium jets over a range of Froude number (~ 0.001 to 
~ I ) w hich is proportional to I /Ri and Rey nolds number 
(~ I to = 100). Strouhal number was correlated as a 
function of inverse Froude number. Oscillations in the 
flow were not observed until a minimum flow rate was 
attained. Measured frequency and the location of the 
vortical structures were dependent on live helium jet exit 
velocity. 

Several experiments by Cetegen on axisymmetric 
helium plumes yielded global oscillatory behavior for a 
range of tube diameters (3.6cm < d < 20cm). The 
velocity measurements in a buoyant plume revealed 
strong buoyant acceleration along the centerline 
followed by a deceleration in live region of toroidal 
vortex fomiation. Cetegen 1 " found live effects of 
external forcing on naturally stable axisymmetric 
buoyant helium-air mixtures. The size of the vortex 
formed was dependent on the forcing frequency in these 
buoyancy-dominated jets < Ri > I). 
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Most recently Yep ci al.' used lire microgravity 
environment of 2.2 seconds drop tower at the NASA 
John H. Cdenn Research Center to find that the flow 
oscillations disappeared in microgravity, demonstrating 
that buoyancy was the origin of self-excited oscillations 
in Earth gravity. Yep et al* provided quantitative 
details of buoyancy-induced changes in the helium 
concentration field at Ri-1.44 using tire Rainbow 
Schlicrcn Deflectomctry (RSDl technique. 

Although the instability in low-density jets is caused by 
buoyancy for RiX).5, in momentum dominated jets, 
however, the origin of the instability still remains 
unknown. Generally, buoyancy effect in momentum 
dominated jets is considered negligible and thus 
ignored. Monkewitz et al.* reported the spectacular 
spreading of unstable treated air jets below a critical 
value of density ratio of approximately 0.73. using 
conventional knife-edge schlicrcn technique. The 
spreading of tire jet was characterized by the formation 
of vortex-ring structures and side-jets. The instability 
was attributed to the vortex-ring formation. 

Kyle et al.® studied the instability and the subsequent 
breakdown of axisytnmetric jets of helium'air mixtures 
diverging into ambient air. Tire oscillating behavior was 
shown to depend upon density ratio, jet diameter and jet 
momentum thickness, and it was independent of the jet 
Rey nolds number within the covered range 1 500 < Re < 
12500. They also noted that tire flow structure repeated 
itself with extreme regularity. The high degree of 
repeatability of the oscillating mode, with strong voitex 
interaction, led to large centerline velocity fluctuations. 
Even though the oscillating mode had high degree of 
repeatability, it did not exhibit large spatial growth rate. 
The instability of the jet was enhanced when the density 
ratio was decreased. 

Richards el al." found helium jets in air to display a 
global or self-excited behavior characterized by intense 
vortex interaction. Flow visualization was perfomred 
using Mie scattering. Concentration measurements at 
several radial locations obtained by an aspirating probe 
indicated the presence of side jet behavior with 
substantially increased level of mixing. Global 
oscillations were repotted for the range of Reynolds 
numbers between 350 and 1400. 

In the present work, we apply the RSD technique to 
measure concentration in sclf-cxeited momentum 
dominated helium jets and obtain quantitative data on 
spatial and temporal evolutions of the instability, and 
length and lime scales of the oscillating mode. 

EXPERIMENTAL APPROAC H 

A compact test rig primarily designed and constructed 
to conduct experiments in microgravity conditions was 
used in this study. An external flow system was coupled 


with the test rig to eliminate limitations regarding the 
maximum flow rate. Flow system was comprised of a 
helium gas cylinder, a pressure regulator, a needle 
valve, a mass flow meter, plastic flexible hose and a 
straight jet tube oriented vertically. 



Figure I. Optical layout of RSD. A) light source 
aperture. B) collimating lens. C) test section. D) 
decollimating lens. E) flat surface mirrors. F) 
rainbow filter. G) high speed camera (All 
dimensions are in millimeters). 

The optical layout of the rainbow schlicrcn apparatus is 
show n in Figure I. It consists of a light source aperture 
w ith 3mm high and 100 pm w idth connected to a 1 50- 
W halogen light source, modified by attaching a 
focusing optics assembly to the exit route of the light, 
through fiber optic cable. 80 mm diameter 310 mm 
focal length collimating achromatic lens. 80 mm 
diameter 1000 mm focal length decollimating 
achromatic lens, two aluminum coated flat surface 
reflecting mirrors. 3 mm wide rainbow filter positioned 
at the focal length of the decollimating lens, and a high 
speed digital camera capable of acquiring images at 
rates up to 2000 frames per second and storing digitally 
in TIFF format at pixel resolution of 384x512. 
Exposure time of 209 psec was used to acquire the 
images. 

In tire rainbow sehlieren apparatus, the transverse ray 
displacement at the filter plane is found from 
measurement of hue in the sehlieren image together 
with a filter calibration curve 11 ' 1 *. The transverse ray 
displacement is related to the angular deflection of a 
light ray by an axisymmetric refractive index field 
given by the follow ing relationship. 



Equation I 


The refractive index difference is found using Abel 
inversion of Equation I given in discrete fomi as 1 ’: 

M 
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For a mixture of tpscs. the refractive index difference U 
given s* 11 

8 = -iLV K XM Equation J 

R-T- 1 

Equation 3 is used to consiiuct a plot between refractive 
index difference and helium mole fraction, assuming 
standard air composing nitrogen to oxygen mole ratio 
of 3.76. 

The uncertainties in the experiments are caused by hue 
measurement uncertainties and instrumental bias eiror. 
The variation of hue values in filter calibration led to 
uncertainties in calculations of deflection angle, 
refractive index and helium concentration. Maximum 
uncertainty in the deficetion angle occurred near the jet 
centerline where the density gradients are small. 
Calculated unceitainty for deflection angle was 5x1 O'* 
degree. The unceitainty in the concentration of helium 
was found to be 8 %. Unceitainty in the jet exit velocity 
was 10.013 ni'S corresponding to uncertainties in 
Reynolds number and Richardson iiumhet of r2 and 
±10"'. respectively. 

RESULTS AND DISC 1 SMON 

In this study, helium concentration measurements are 
presented for the follow ing two test cases w ith d= I ‘>.05 
mm. 


rest 

Re 

Ri 

1 

91(1 

0.034 

2 

1250 

0 018 


First, w e w ill discuss oscillation characteiisties of self- 
evened momentum dominated jets. Power spectra taken 
at various locations in the flow field will be presented 
in this pan. Then, we w ill examine the instantaneous 
contours of helium mole fraction duiing an oscillation 
cycle, and show the repetitiveness of the flow and the 
vortical structures. Finally, tlx? temporal evolution plots 
of helium concentrations will be presented and 
discuswd. 

Owillalinc C luiacteri>!ics and Global Feature-' 

Sample images for the two test eases are shown in 
Figure 2. The hue in the background and jet center is 
caused by the undofleetcd light rays where die density 
gradients are small or absent. Hue othei than the 
background hue is caused by light ray deflection due to 
density gradients in the medium. Power spectra at 
vaiious locations in the near flow field were generated 
using a Fast Fourier transform (FFT) algorithm 



a n 


Figure 2. Rainbow vchliemi Images. A) test case I. 
B) test case 2. 

To obtain flic power spectra, first the time-space images 
were generated at desired axial locations using 20-18 
consecutive rainbow sehlieren images. A typical lime- 
space image for test ease I is shown in Figure 3. The 
veitical direction in Figure 3 indicates the radial 
coordinate and the horizontal direction represents the 
time. For improved visualization, only 0.5 sec time 
interval is sltown. After generating die space-time, the 
hue values at desired radial cooulinates were used as 
the input for FFT analy sis to obtain the power spectra. 



Figure 3. Time-space image at z/d-l for lest case I. 

Figure 4 shows power spectra at r/d=0.45 for different 
axial locations (zd=0.05. 0.5 and 1.0) for test case I. 
Results depict a dominant oscillation frequency of 125 
Hz at all axial locations. Figure 4 shows that the 
oscillation amplitude increased in the flow direction. 
The increase in die oscillation amplitude is attributed to 
the growth of the vortical structure affecting die flow 
field. Figure 5 shows power spectra at z'd- 1 .0 for 
dilTcrent radial locations (i d =0.55. 0.6 and 0.(>5). 
Again, a dominant oscillation frequency of 125 11/ was 
obtained at all radial locations signifying die oscillation 
mode is global and independent of spatial coordinates. 
Hie highest oscillation amplitude was obtained at 
r/d-0.6. 


3 

American Instituie of Aeronautics and Astronautics 





Figure 5. Power spectra at u'd-l for test case I. 


Figure 6 shows power spectra for test case 2 at z/d=I.O 
and r/d-0.4 to demonstrate the oscillatory behavior. A 
global oscillation frequency of ISO Hz was found for 
this case. 

Instantaneous Contours of Helium Mole Fraction 

Instantaneous contours of the helium concentration 
were generated to explore the oscillatory behavior of 



Frequency mi) 


Figure 6. Power spectral at z/d=l and r/d=0.4 for 
test case 2 . 


the flow during a cycle. Results are shown in Figure 7 
for test case I. Owing to the axisymmetric flow, only 
the one side of the jet is shown. At t=0 ms and t-8 ms 
(not shown), the contours are identical to each other 
suggesting the oscillation period of 8 ms or oscillation 
frequency of 125 Hz. In Figure 7. the concentration 
field is changing throughout the cycle to depict the 
instability of the flow. The instability is caused by the 
vortical structure which evolves and convecls 
downstream during the oscillation cycle. The vortex 
initiation and its propagation to the downstream can be 
examined from Figure 7. At t=0 ms, a vortex forms 
around z'd=0.25 at>d r'd=0.48 as seen from the 
indentation on tire 70 % concentration level. At t- 1 ms. 
the vortex located at z/d-0.40 has expanded in both 
axial and radial directions. At t=2ms. the vortex is 
located around z'd=0.S5. Subsequently, as the vortex 
convecls dow nstream it affects a greater portion of the 
jet. Figure 8 shows the instantaneous contours of 
helium concentrations for test case 2. Again, the overall 
trends are similar to those observed at a lower Reynolds 
number (test case 1 ). 
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Temooial Evolution of Helium Mule Fraction 

In ihc previous section, the oscillatory behavior of the 
flow was discussed using instantaneous helium 
concentration contours for one oscillation cycle. 
Temporal evolution plots of helium concentration 
prov ide detailed information about the flow structure by 
containing multiple cycles. Figure 9 shows the temporal 
evolution plots of helium concentration for the test case 
I. 








Figure 9. Temporal evolution of helium mole 
fraction at different axial locations for test case I 
(left side) and for test ease 2 (right side). 


Time traces were shown for 0.05 sec interval for one 
side of the jet centerline. Changes in the Itclium 
concentrations at different axial and radial locations can 
be monitored from plots in Figure 9. Results show 
highly repetitive helium concentration at all locations 
indicating a global oscillation mode. The oscillation 
frequency based on time interval between consecutive 
peaks was 125 Hz. which agrees with the value 
obtained using FFT analysis. Tire oscillation amplitude 
increases in the flow direction confirming the results 
obtained from power spectra plots. Small amplitude 
oscillations were observed near the tube exit at 
zd-0.05. At zd-0.5. the oscillations are dominant as 
seen from increased concentration amplitudes and 
helium has diffused up to r/d-0.71. At z'd= I . the 
concentration amplitude has increased about 54 % 
compared to that at zd-0.5. At z'd-1.5. the amplitude 
of oscillation is very high. The oscillations penetrate 
deeper into jet core. Helium has diffused further into air 
as seen by 10 % concentration level located at r'd=0.81. 
In Figure 9. the temporal evolution of helium 
concentrations are show n for the test case 2. Results are 
similar to those obtained for test case I. 

CONCLUSIONS 

Flow structure of self-excited, laminar, axisymmctric, 
momentum-dominated helium jets discharged vertically 
into ambient air was investigated for two test cases 
using high-speed rainbow schlicren deflectometry 
technique. Power spectra plots revealed intense 
oscillations marked by a dominant frequency 
independent of the radial and axial coordinates. The 
increase in peak power with flow direction was also 
observed from these plots. Instantaneous contours of 
helium mole percentage show ed initiation of the v ortex 
and its interaction with the flow field duiing an 
oscillation cycle. Temporal evolution plots of helium 
concentrations provided information on variations in 
the jet structure with flow direction. 
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NOMENCLATURE 

d Tube Inside diameter 

f Frequency! Hz) 

g Gravitational Acceleration 

r Radial Distance from Centerline 

y Projected Location 

D Geometric Coefficients 

M Species Molecular Weight 

P Atmospheric Pressure (N/m 5 ) 
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R Universal Gas Constant 

Re Reynolds Number 

Ri Richaulson Number 

5 Density Ratio 

T Temperature 

U Mean Jet Exit Velocity (nVs> 

X Species Mole Fraction 

Greek Symbols 

£ Deflection Any let radian ) 

»; Dalc-Gladstone Constant 

6 Refractive Index Difference 

t| Refractive Index Normalized by that of 

Surroundings 
p Densityfkg/ m 1 ) 

0 Momentum Thickness 

v Kinematic V'iscosity (nr/s) 

Subscripts 

<*> Ambient Air 

1 Projected Radial Location 

j Species 
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